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We investigate the presence of defects in systems described by real scalar field in (D,l) spacetime 
dimensions. We show that when the potential assumes specific form, there are models which support 
stable global defects for D arbitrary. We also show how to find first-order differential equations that 
solve the equations of motion, and how to solve models in D dimensions via soluble problems in 
D=l. We illustrate the procedure examining specific models and finding explicit solutions. 

PACS numbers: ll.10.Lm, 11.27.+d, 98.80.Cq 
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The search for defect structures of topological nature 
is of direct interest to high energy physics, in particular 
to gravity in warped spacetimes involving D spatial extra 
dimensions. Very recently, a great deal of attention has 
been given to scalar fields coupled to gravity in (4, 1) 
dimensions Q, 0, El 0, 0- Our interest here is related 
to Ref. 0, which deals with critical behavior of thick 
branes induced at high temperature, and to Refs. 
which study the coupling of scalar and other fields to 
gravity in warped spacetimes involving two or more extra 
dimensions. 

These specific investigations have motivated us to 
study defect solutions in models involving scalar field in 
(D, 1) spacetime dimensions. To do this, however, we 
have to circumvent a theorem 0, 0, 0] , which states 
that models described by a single real scalar field can- 
not support topological defects, unless we work in (1,1) 
space-time dimensions. To evade this problem, in the 
present letter we consider models described by the La- 
grange density C = (1/2)9^09^0 — U{x 2 ]<f)), where the 
potential is U(x 2 ;4>) — f(x 2 )V((t>) 1 x 2 = x^x^, and (f> 
is a real scalar field. The metric is (+,—,...,—), with 
= (xq, X\, x%i Xd) in (-0,1) space-time dimensions. 
V(4>) has the form V{<j>) = W|/2, where W = W((j>) is 
a smooth function of 0, and W$ — dW/d<fr. We sup- 
pose that is a critical point of V, such that V{<j>) = 0. 
This generalization is different from the extensions one 
usually considers to evade the aforementined problem, 
which include for instance constraints in the scalar fields 
and/or the presence of fields with nonzero spin - see, 
e.g Ref. [ljj, and other specific works on the subject 
|l4Lll5| . Potentials of the above form appear for instance 
in the Gross-Pitaevski equation, which finds applications 
in several branches of physics - see, e.g., Ref. [16|. Other 
recent examples in (1,1) dimensions include Ref. [l7j . 
which deals with the dynamics of embedded kinks, and 
Refs. ,18|, which describe scalar field in distinct back- 
grounds. 

In higher dimensions, the factor l/r N that we intro- 
duce in Eq. gives rise to an effective model, which 
comes from a more fundamental theory. To make this 
point clear, we consider the model C\ = d^d^cj) — 
f({p)F flv F t11 ', which is a simplified Abelian version of the 



color dielectric model 0] in the absence of fermions - 
see, e.g., Ref. JjjjJ. This model describes coupling be- 
tween the real scalar field and the gauge field , through 
the dielectric function /(</>). Here F^ v = d ii A v — d v A^ 
is the gauge field strenght. This model shows that for 
spherically symmetric static configurations in the elec- 
tric sector, the equation of motion for the matter field 
is V0 + {df/d^E 2 = 0, where E = (E r , 0, 0) is the 
electric field. The use of the equation of motion for the 
electric field provides an effective description of the mat- 
ter field <j), in which the equation of motion exactly repro- 
duces |2l| the equation of motion that we will be inves- 
tigating below, under the choice /(</>) = 1/V((f>). As one 
knows, color dielectric models may provide effective de- 
scriptions for non-perturbative QCD - see, e.g., Ref. 
for a very recent investigation and for related works. 

We now focus attention on the effective model. We 
investigate the presence of static solutions considering 
the specific potential 
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where r = (xf + x\ + ■ ■ ■ + x 2 D ) 1/2 . We write 4> = <f>{r) to 
obtain the equation of motion V 2 = (l/2r N )W ! f,W^. 
The energy density of the field configuration has the form 
£ = (1/2) (V0) 2 + (l/2r N )W$. We suppose that (j> = 
4>{t) solves the equation of motion. It has total energy 
E D = E® + Ep, which splits into gradient and potential 
portions. We make the change <f>(r) — » <fi x (r) — 4>{\r) to 
get to the conditions 
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and (2-D){l-D)Ef + {N -D){N -D-l)E% > 0, 
to make the field configuration stable. Since E£ > and 
Ep > 0, these conditions impose restrictions on both N 
and D. 

An important case is D = 1, and here the value N = 
makes Eg = Ei, that is, the energy is equally shared be- 
tween gradient andpotential portions. This is the stan- 
dard result: from [lOl llll Il2| we see that for the special 
case N = 0, there exist stable solutions only in D = 1. 
However, the form of the potential Q is peculiar, and it 
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allows obtaining several other cases which support defect 
solutions. In two spatial dimensions, for D = 2 one gets 
N = 2, but this gives no further relation between the 
gradient and potential portions of the energy. For D > 3 
we get that N = 2D - 2 to make Ef = . 

We investigate t he p ossibility of obtaining a Bogo- 
mol'nyi bound [13, l24j for the energy of static config- 
urations in the present investigation. The case D = 1 
is standard, so we deal with D > 2. We suppose that 
the static solutions engender spherical symmetry, that 
is, we consider tf> = <j)(r), with r G [0, oo). The en- 
ergy can be written as E D > Cljj\AW\, where AW = 
±W[{<j>{r -> oo)}] =F W[{(j){r = 0)}] and Q D is the D- 
dependent angular factor. This result is obtained if and 
only if N = 2D — 2, and the energy is minimized to 
E D = fl£)\AW\ for static and radial field configurations 
which solve the first order equations 
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This is the Bogomol'nyi bound, now generalized to 
models for scalar field that live in D spatial dimen- 
sions. The value N = 2D — 2 leads to defect solu- 
tions of the BPS type, which obey first order equations 
and have energy evenly split into gradient and poten- 
tial portions. For D = 2 one gets N = 2, and in 
this case the model engenders scale invariance. We can 
show that solutions of the above first order equations 
solve the equation of motion JU) for potentials given by 
Eq. QJ. Also, we follow Ref. |25j and introduce the ra- 
tio R = r D ~ 1 (d(j)/dr)/W c f > . For field configuration that 
obey 4>(0) — 4> an d lim r ^o dcj>/dr — > one can show that 
solutions of the equation of motion also solve the first or- 
der equations © ■ This extends the result of Ref. [25j to 
the present investigation: it shows that the equation of 
motion completely factorizes into the two first order 
equations Q. 

The equation of motion for <f> = 4>{r) is 
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The first order equation is given by (J3J, and their solu- 
tions solve the equation of motion and are stable against 
radial, time-dependent fluctuations. To see this we con- 
sider cj)(r, t) = 4>{r) + J^k Vk{r) cos(wk t). For small fluc- 
tuations we get Hr/k 
be written as 



tu|?7fe, where the Hamiltonian can 



H 



1 



„2-D-2 



(5) 

It is non-negative, and the lowest bound state is the zero 
mode, which obeys r D ~ 1 dr] /dr = ±W^<prio- This gives 
7 ?o( r ) = c exp (± J drr Y ~ D W^) , where c is the normal- 
ization constant, which usually exists only for one of the 



two sign possibilities. We can also write 770 (r) = cW^,, 
which is another way to write the zero mode. 

We investigate the presence of defect structures turn- 
ing attention to specific models in D = 1, D = 2, and 
D > 3. We first consider the case D = 1. We choose 
N = 0. The equation of motion is (d 2 (j>/dx 2 ) = W^. 
We are searching for solutions that obey the boundary 
conditions lim^^-oo <fi(x) — > <f> and liTa x ^ 00 (d(/)/dx) — > 
0, where is a critical point of the potential, obey- 
ing V{4>) = 0. In this case one can write (d<j)/dx) — 
±(dW/d(f>); see [isf. In D = 1 one usually introduces the 
conserved current j M = e^d u <j). We see that p — d<p/dx 
and so gives the energy density of the field configu- 
ration (26| . Thus, we introduce Qt — dxp 2 as the 
topological charge, which is exactly the total energy of 
the solution. 

We exemplify the case D = 1 with the family of models 
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The parameter p is real, and it is related to the way the 
field self- interacts. These models are well-defined for p 
odd, p = 1,3,5... For p = 1 we get to the standard 4> 4 
theory. For p = 3, 5, ... we have new models, presenting 
potentials which support minima at <fi = and ±1. For p 
odd the classical bosonic masses at the asymmetric min- 
ima 4> = ±1 are given by m 2 = 4/p 2 . For p — 3,5,... 
another minimum appear at <p = 0. However, the classi- 
cal mass at this symmetric minimum diverges, signalling 
that 4> = does not define a true perturbative ground 
state for the system. 

We consider p odd. The first-order equations are 
d(f)/dx — ±4>~ =F0~, which have solutions <j>± '(x) = 
± tanh p (x/p). We consider the center of the defect at 
5 = 0, for simplicity. Their energies are given by 
e[ p I = Ap/ (4p 2 — 1) , and we plot some of them in Fig. [1] . 
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FIG. 1: Defect solutions for models with p — 3 and 5, respec- 
tively full and dashed-dotted lines. 

We see that solutions for p = 3, 5, connect the min- 
ima <j> = ±1, passing through the symmetric minimum 
at <j) — with vanishing derivative. They are new struc- 
tures, which solve first-order equations, and we call them 
2-kink defects since they seem to be composed of two 
standard kinks, symmetrically separated by a distance 
which is proportional to p, the parameter that specifies 
the potential. To see this, we notice that the zero modes 
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are given by ?7q 



(i,p) 



c^-rt tanh p-1 (x/jj) sech 2 (x/p), 
where c^ 1 '^ = [(4p 2 — l)/^] 1 / 2 - we use rjQ D ' p ^ to rep- 
resent the zero modes. These zero modes concentrate 
around two symmetric points, which identify each one 
of the two standard kinks. Defect structures similar to 
the above 2-kink defects have been studied in the recent 
past, for instance in Ref. (2^, for solutions of the equa- 
tions of motion in a Z(3)-symmctric model, and also in 
Ref. |28( , in a supersymmetric theory, for solutions which 
solve first-order equations. We have studied the V p {4>) 
model coupled to gravity in warped spacetimes, in (4, 1) 
dimensions [2^|. The results show that it also induces 
critical behavior similar to Ref. Q, but now the driving 
parameter is p, which indicates the way the scalar field 
self-interacts, and not the temperature anymore. 

The cases with p even are different. The case p = 2 is 
special: it gives the potential V((p) = cp/2 — <p 2 + 4> 3 /2, 
which supports the nontopological or lumplike solution 
<j>l' 2 '(x) — tanh 2 (a;/2). The lumplike solution is unsta- 
ble, as we can see from the zero mode, which is pro- 
portional to tanh(x/2)sech 2 (x/2): the zero mode has a 
node, so there must be a lower (negative energy) eigen- 
value. In the present case, the tachyonic eigenvalue can 
be calculated exactly since the quantum-mechanical po- 
tential associated with stability of the lumplike configura- 
tion has the form U(x) = 1 — 3 sech 2 (a;/2). This potential 
supports three bound states, the first being a tachyonic 
eigenfunction with eigenvalue Wq = —5/4, the second the 
zero mode, w 2 — 0, and the third a positive energy bound 
state with w 2 = 3/4. As one knows, tachyons appear in 
String Theory [3(], E01 an d this has br oug ht renewed in- 
terest on the subject - see, e.g., Refs. Ullillli- 

The other cases for p even are p = 4, 6, ... These cases 
require that > in Eq. 10, but we can also change 
6 — » — <t> in Eq. and consider <fi < 0. We investigate 
the case 4> £ [0, oo); reflection symmetry leads to the 
other case. We notice that the origin is also a minimum, 
with null derivative. These models also support topolog- 
ical defects, in the form cf>^-' p '(x) — tanh p (a;/p) (x > 0), 
with energies e[ p I = 2p/{Ap 2 — 1) for p — 4, 6, ... These 
solutions solve the first order equation dcj>/dx — W^,; 
the other equation dip/dx — — is solved by <p{x) — 
-tanh p (a;/p) (x < 0). 

We now consider D — 2, and N — 2. The equation 
of motion is V 2 = (l/r 2 )W^W^0. We search for solu- 
tion <p(r), which only depends on the radial coordinate, 
obeying 4>(0) = 4> and lim r ^o d<j)/dr — > 0. In this case we 
get 
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We write dx = ±r~ 1 dr to get d 2 (pjdx 2 = W^W^. This 
result maps the D — 2 model into the D = 1 model. 

We see that r = exp(±.x), which shows that the full 
line x £ (—00,00) is mapped to r £ [0, 00). We notice 



that since the center of the defect is arbitrary, so is the 
point r = 1 in D = 2; thus the solution introduces no 
fundamental scale, in accordance with the scale symme- 
try that the model engenders. If one uses the model © 
with p odd to define the potential in this case, we get the 
solutions 



Their energies are = 8irp/(4p 2 — 1). In Fig. [2] we 
depict the defect solution for p — 1. The corresponding 
zero mode is given by % ( r ) = V '^2 / '^[r 2 / \r 2 + l) 2 ]. 
This zero mode binds around the circle where the defect 
solution vanishes, as we show in Fig. [2], where we also 
plot p£' 1] (r) =32r 4 /7r(r 2 + l) 4 . 
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FIG. 2: Radial defect solutions for D — 2, in the case p = 1. 
The dashed-dotted line represents the density Pq 2 ' 1 ' (r) of the 
corresponding zero mode. 

The last case is D > 3, with N = 2D — 2. The equation 
of motion for <h = <b(r) is 
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dr 
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We write dx = ±r 1 ~ D dr to get d 2 4>/dx 2 = W^W^ and 
again, we map the model into a one-dimensional problem. 
We solve dx = ±r 1 - D dr to get x = ^A 2 - ^ /{D - 2). 
This shows that x is now in (—00, 0] or [0, 00), and so we 
have to use the upper sign for x < 0, or the lower sign for 
x > 0. We can use the model of Eq. © for p = 4, 6, ... to 
solve the D-dimensional problem with ./V = 2D — 2. In 
this case the solutions are 



(r) = tanh J 
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for D = 3,4,... Their energies are given by E£ 
£Id[2p/(4:P 2 — 1)], and in Fig. [3] we depict the solution 
for D = 3. 

In the case of D — 3 and p = 4 the zero mode is given 
by Vo 3 ' 4 \ r ) = c(M) tanh 3 (l/4r)sech 2 (l/4r); we could not 
find the normalization factor explicitly in this case. We 
plot p^' 4 \r) in Fig. [3] to show that the zero mode 
binds at the skin of the defect solution, concentrating 
around the radius R of the defect, which is given by 
R = l/4arctanh[(l/2) 1 / 4 ]. 
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The solutions that we have found have a central core, 
and a skin which depends on the parameters that spec- 
ify the potential of the model. They are stable, well 
distinct of other known defects such as the bubbles 
formed from unstable domain walls - see for instance 
Refs. [35l US 113, HU • Also, they are neutral structures, 
and may contribute to curve spacetime, and to affect 
cosmic evolution. They may become charged if charged 
bosons and/or fermions bind to them. For instance, if 
one couples Dirac fermions with the Yukawa coupling 
Y{(j>) = r 1 "' Wtptf,, the fermionic zero modes are similar 
to the bosonic zero modes that we have just presented. 
In Ref. [1H we investigate these and other issues in detail. 




FIG. 3: Radial defect solutions for D — 3, in the case p = 4. 
The dashed-dotted line shows (1/10 of) the density of the 
corresponding zero mode. 

We end this letter recalling that we have solved the 



equation of motion using 



,1-D 



dr = ±dx. Although 



this identification works very naturally for the first or- 
der equations, we used the equation of motion to present 
a more general investigation, which may be extended to 
tachyons in D spatial dimensions. For instance, from the 
lumplike solution for D — 1 and p — 2 we can obtain 
another lumplike solution <^>| 2 ' 2 ^(r) = (r — l) 2 /( r + I) 2 ; 
which is valid for D = 2. Another issue concerns the 
identification of the topological behavior of the above de- 
fect structures. We do this introducing the (generalized 
current-like) tensor jw-vo = £ ^iM2---mdmd+i^ d+i 0. 
It obeys d IH j ttltl3 " ,tlD = for i = 1, 2, • • • , D, which 
means that the quantities p 1112 '" 10 - 1 = j°nvs— con _ 
stitute a family of D distinct conserved (generalized 
charge) densities. We introduce the scalar quantity 
pl = Plll2 ... lD _ lP l ^- l °-i = (~l) D (D - l)\{dcj>/drf, 
which generalizes the standard result, obtained with 
jfi _ £ ^d v< j) - see the reasoning above Eq. ©. Thus, 
we define the topological charge as — J drp 2 D = 
(— 1) D (D — ly.flD^W, which exposes the topological be- 
havior of the new global defect structures that we have 
just found. 

We would like to thank R. Jackiw for drawing our at- 
tention to Ref. [13 • We also thank CAPES, CNPq, PRO- 
CAD and PRONEX for partial support. 
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